Substitutional doping of Cu in diamond: Mott physics with p orbitals 
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Discovery of superconductivity in the impurity band formed by heavy doping of boron into diamond (C:B) as 
well as doping of boron into silicon (Si:B) has provided a rout for the possibility of new families of supercon- 
ducting materials. Motivated by the special role played by copper atoms in high temperature superconducting 
materials where essentially Cu d orbitals are responsible for a variety of congelation induced phases, in this paper 
we investigate the effect of substitutional doping of Cu into diamond. Our extensive first principle calculations 
averaged over various geometries based on density functional theory, indicates the formation of a mid-gap band, 
which mainly arises from the t2g and 4p states of Cu. For impurity concentrations of more than ~ 1%, the 
effect of 2p bands of neighboring carbon atoms can be ignored. Based on our detailed analysis, we suggest a 
two band model for the mid-gap states consisting of a quarter-filled hole like t2g band, and a half-filled band of 
4p states. Increasing the concentration of the Cu impurity beyond ~ 5%, completely closes the spectral gap of 
the host diamond. 

PACS numbers: 71.55.-i, 71.23.-k, 78.30.Am 
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INTRODUCTION 

Diamond is the best thermal conductor |[ll|2], and the hard- 
est known material with semiconducting properties. Both 
pristine form as well as doped films of diamond present 
unique ground for high-temperature and high-power applica- 
tions in the industries Recent discovery of supercon- 
ductivity in heavily boron doped diamond, has revived the in- 
terest in impurity band formation from a fundamental point 
of view Jsl Bl Similar system of boron doped silicon was 
also studied ITHSt]. Synthesis of boron doped diamond was 
first achieved by Ekimov and coworkers jsd and subsequently 
by Takano et. at. |@]. They used chemical vapor deposition 
(CVD) for the synthesis of B-doped diamond. At low doping 
concentrations of typically 10^^ — 10^^ cm^'^, boron atoms in- 
ject acceptor levels in the electronic energy spectrum of host 
diamond, hence producing a p type semiconductor ifioll . For 
heavier dopings ~^[B]/[C]> 5000 ppm in the gas phase, one 
achieves metallic conductivity in diamond Eventually 
increasing the doping beyond n > 10^^ cm^'^, i.e. few per- 
cents, it superconducts at low temperatures [5, 6]. For higher 
concentrations, the dopant atoms come closer to each other 
and start to overlap. Therefore the imp urity levels are broad- 
ened into authentic impurity bands 112111311 . which are respon- 
sible for metallic and superconducting properties Is. 13- 181. 

Various forms of defects play crucial role in the physical 
properties of diamond. For example the fascinating colors 
of diamond as gem stone, is due to the ~^ppm concentration 
of vacancies fll9tl . Such vacancies at much higher concentra- 
tion of the order of few percents, can give rise to an impurity 
bands inside the diamond gap, which can provide larger num- 
ber of states than some of the elements from the third and 
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fifth column of the periodic table of elements 112011 . Such mid- 
gap bands of the vacant diamond can be attributed to a local- 
ized Wannier wave function composed from the 2p states of 
neighboring carbon atoms. Therefore, previous studies of the 
possible conductive impurity bands in diamond and/or silicon 
hosts, have been mainly concerned with the substitution of a 
p-type orbital for the host p-states. 

In this paper, we are interested in the role of substitution by 
a 3d element. We choose to study the effect of copper sub- 
stitution in the diamond host. Why is Cu interesting from the 
impurity band formation point of view? First of all, the same 
high-temperature, high-pressure techniques employed for the 
synthesis of boron doped diamond can also be applied to at- 
tempt a substitution of Cu atoms into the diamond host 12 ih . 
From theoretical point of view, although Cu is a 3d transition 
metal, it does not have magnetic complications, and hence one 
can avoid possible magnetic ground states at low temperatures 
which tend to compete against superconducting correlations. 
Choice of Cu also avoids the possibility of local moment for- 
mation. On the other hand 3d^ "atomic" configuration of Cu 
is hoped to provide a hole which plays an important role in 
the high temperature superconducting (HTSC) materials 12211 . 
The nominal 3d hole of atomic Cu in HTSC systems is shared 
between neighboring oxygen 2p and Cu 3d states, giving rise 
to the so called Zhang-Rice singlet state 112 311 which caiTies the 
hole. 

In this paper, we study the formation of impurity bands in 
Cu doped diamond in detail and we find that 4p states of Cu 
along with its t2g orbitals are responsible for the formation of 
impurity band at Cu concentrations 1 — 5%. At lower concen- 
trations 2p orbitals of neighboring carbon atoms also partic- 
ipate in the formation of mid-gap metallic band. Depending 
on the experimental conditions, either a combination of most 
probable geometries is realized, or by slow cooling the system 
will be given enough time to perform a self-averaging over the 
phase space of possible geometries. Therefore we undertake 
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extensive averaging over distinct geometries. The main role of 
averaging is to restore the symmetries broken by the specific 
realization, e.g. most probable geometry (MPG). For com- 
parison, we will provide MPG and geometry averaged (OA) 
results together 



METHOD OF CALCULATION 

The calculations of electronic structure are performed with 
the plane-wave pseudopotential 12411 of DFT using the PBE 
functional 125] for exchange and correlation as implemented 
in the QUANTUM-ESPRESSO package. Ultra-soft pseu- 
dopotentials are used with a 30 Ry (230 Ry) cutoff for the 
expansion of the wave-functions (charge density). We have 
constructed 2x2x2 supercell containing 64 atoms with 1, 2 
and 3 Cu impurity, leading to doping concentrations of 1.6%, 
3.1% and 4.7%, respectively. The sampling of Brillouin zone 
(BZ) for structural relaxation and DOS calculations are per- 
formed on 2 X 2 X 2 and 4x4x4 grid of fc-points, respectively. 
We first minimized the energy of undoped diamond with re- 
spect to cell size to calculate the lattice constant. The obtained 
value was very close to the experimental result. All structures 
were relaxed by BEGS algorithm ll26ll . 

Two types of doping can be considered: substitutional and 
interstitial. In the case of B-doped diamond, interstitial dop- 
ing is unstable ||27il . Since our study is motivated by B- 
doped diamond, we only consider substitutional doping of 
copper. We randomly substituted Cu impurity into the dia- 
mond structure. Since there are many different ways of in- 
troducing impurity into the host carbon lattice, one needs to 
average over sufficient number of realizations of randomness. 
With m impurities in a host of n lattice points, there are 
C{n, m) = (n~my.m,\ different ways of substitution. With 
n = 64 carbon atoms considered here, the most trivial case 
corresponds to m = 1 impurity. By translational symmetry 
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(periodic boundary conditions), it does not matter where we 
place the single impurity. Therefore the average DOS cor- 
responding to 1 impurity, is identical to any of the random 
realizations. Next in the case of m = 2 impurities, there are 
in principle C(64, 2) = 2016 ways of placing the Cu impuri- 
ties. However, since the relative distance between the impuri- 
ties is the physically discriminating factor, many of the above 
configurations will be essentially identical by translation and 
rotation symmetries. Therefore the electronic structure calcu- 
lations must be run for geometrically different configurations. 
We have used a computer program to count the number of ge- 
ometrically independent configurations which in the case of 
m = 2 impurities gives 23. Therefore one has to perform the 
electronic structure for all 23 distinct geometries, taking into 
account the multiplicity Wg of each geometry g. In the case 
of m = 3 impurities, we have C(64, 3) ~ 41664 nominally 
different configurations. Again our counting algorithm gives 
430 distinct geometries, where each geometry is meant to rep- 
resents a distinct set of 3 mutual distances between the impu- 
rity atoms. We used 60 most probable geometries to calculate 
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FIG. 1: (Color online) Total density of states for three considered 
concentrations contrasted against pure diamond in background. Left 
and right panels correspond to geometry averaged and most probable 
geometries, respectively. 



the averages. Note that the above 60 geometries cover more 
than half of the total 41664 choices. We have checked that, 
the average DOS obtained for 60 most probable geometries 
can be more or less achieved even with 10 most likely geome- 
tries. The geometry averaged values of various quantities X 
such as p(a;), Ep can be calculated as follows; 



^GA = ^WgXg, 



(1) 



where g runs over geometries, and Wg is its relative weight. 



METALIZATION OF DIAMOND 

In Eig.[T|we have plotted the DOS of doped diamond lattice 
with (left) and without (right) averaging over geometries. The 
doping rates in the figure correspond to 1, 2, 3 Cu impurities 
in a host of 64 carbon atoms, from bottom to top, respectively. 
The presence of Cu impurities, affects the entire spectrum (not 
shown here). In this figure we have focused on the region of 
spectrum corresponding to gap in the pristine diamond. By 
increasing the doping rate from 1.6 to 4.7 percent, the devia- 
tion in the DOS of Cu-substituted diamond from the DOS of 
perfect diamond becomes more and more significant. 

The most spectacular effect of substituting Cu in the in- 
sulating host of diamond, is the formation of mid-gap impu- 
rity band. In the right panel of Eig.[Tl DOS corresponding to 
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FIG. 2: (Color online) LDOS for s, p and d orbitals of Cu corre- 
sponding to three concentrations considered. Fermi level Ef is in- 
dicated by vertical dashed line. MPG values on the right panel are 
given for comparison. As can be seen the contribution of 4s orbitals 
to the impurity band is negligible. 
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FIG. 3: (Color online) Copper d-orbital splitting for three considered 
concentrations in diamond crystalline field. The major contribution 
to the impurity band comes from t2g subspace of d-orbitals. 



MPG is shown. By averaging over different geometries, the 
corresponding geometry averaged DOS in the left panel will 
be obtained. As can be seen in the figure, for both GA and 
MPG, the width of the impurity band increases with increas- 
ing the Cu concentration. At 4.7% Cu substitution, the gap 
completely disappears. The averaged Fermi energy, Ep is in- 
dicated with vertical dashed line. As it is evident, although 
the most probable geometry in the case of 2/64 impurity con- 
centration gives a little gap in the spectrum, other geometries 
with high probability will fill in that small gap and the result- 
ing GA density of states will be a metallic band. This indicates 
that the best known insulator, namely diamond, can be metal- 
ized by few percent Cu substitution. The question we would 
like to investigate in detail in this paper is, what is the origin 
of the metallic band? 



Role of Cu atoms 

To consider this question, we first examine the occupation 
of various orbitals. We find that unlike the atomic configura- 
tion 4:S^3dP of Cu, in the presence of carbon neighbors, the 
occupancy of 4s orbitals is changed as 4s^ — )• 4s°'^. The 
missing 1.5 electron changes the — >■ Sd^'^, 4p° — >■ 4p^. 
Hence 4p orbitals are pushed down to acquire occupancy of 



^ 1. We find that the above Lowdin charges are more or less 
independent of the geometry, as well as concentration of Cu 
impurities. This can be qualitatively seen in Fig.|2] where lo- 
cal density of states (LDOS) for s, p and d orbitals of Cu are 
plotted. Therefore the relevant orbitals of Cu are both Ap, as 
well as a subspace of 3d states. 

To further identify which subspace of Cu orbitals con- 
tribute to impurity band, in Fig. [3] we resolve the contribution 
of t2g and eg sub-bands. In the case of one Cu atom, the impu- 
rity atom feels a perfect diamond environment; thereby giving 
three-fold degenerate t2g band and two-fold Cg bands. Note 
that Cg bands are located at lower energies with respect to t2g 
bands. By increasing the number of Cu atoms, the above crys- 
talline degeneracies start to be lifted, as the crystalline envi- 
ronment at higher Cu concentrations will deviate from perfect 
diamond symmetry. At the same time the average distance 
between Cu atoms decreases. Especially in the case of three 
Cu atoms, diamond crystalline environment is totally lost in 
the DOS of MPG. Even after averaging over geometries, sub- 
stantial degeneracy lifting in each t2g and Cg sub-bands can be 
observed. The important result from Fig.[3]is that in all three 
considered concentrations, and for both MPG and GA situ- 
ations, only t2g sub-band of the d-orbitals contribute to the 
impurity band, while the eg sub-band will give energy levels 
inside the valence band. 
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FIG. 4: (Color online) Spatially resolved DOS at the Cu site. Right 
panel corresponds to MPG, and the left panel includes averaging over 
geometries. At high concentrations the dominant contribution comes 
from Cu atom itself, while at lower concentration nearest neighbor 
carbon atoms have substantial contribution to the impurity band. 
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FIG. 5: (Color online) Contribution of first and second neighbors 
(NN and NNN) to the DOS of the impurity band. Right panel corre- 
spond to MPG and the left panel includes averaging over geometries. 
As can be seen s orbitals of neighboring carbon atoms do not con- 
tribute to the impurity band. The relative contribution of neighboring 
p orbitals to the Cu atomic orbitals become significant in low concen- 
tration limit. 



Role of neighboring carbon atoms 



Now let us concentrate on the effect of neighboring carbon 
atoms on the formation of the impurity band. In Fig. |4]we 
plot DOS corresponding to three concentrations for MPG and 
GA. In terms of the relative weights of neighbors, both MPG 
and GA give similar results. Note that for one Cu impurity, 
MPG and GA quantities must always be the same. By moving 
from bottom to top (increasing the concentration) in Fig.|4]one 
notices that only at low concentration nearest neighbor (NN) 
carbon atoms contribute significantly. By increasing the con- 
centration, contribution of NN carbon atoms compared to the 
Cu atoms is reduced and becomes as low as the next nearest 
neighbor (NNN) contributions. This can be attributed to the 
fact that at higher Cu concentration, average Cu-Cu distance 
is reduced and the Cu-Cu hopping matrix elements become 
dominant. At lower concentrations in addition to Cu-Cu hop- 
pings, there is also substantial hopping amplitude between the 
dangling orbitals of neighboring carbon atoms. This can be 
thought of hopping between appropriate Wannier orbitals lo- 
calized around Cu atoms which is composed of the dangling 
2p orbitals of neighboring carbon atoms. 

Based on a naive sp^ picture of bonding in perfect diamond, 
one expects 2s and 2p orbitals of neighboring carbon atoms 
to have comparable contribution to the formation of impurity 
bands. However, Cu substitution drastically changes this pic- 



ture. As can be seen in Fig. HI in both MPG and GA situations 
s orbitals are irrelevant and the contribution of NN and NNN 
carbon atoms comes from their 2p states. These contributions 
become comparable to local DOS arising from Cu atoms only 
in low Cu concentration limit. When comparing Figs. S] |5] 
one has to note the scales. Note that 1/64 concentration can 
also be reproduced in the simulation with two Cu atoms in 
host supercell of 128 carbon sites. Therefore the growth in 
the relative weight of neighboring p orbitals can be partly due 
to finite size artifacts of the simulation, which of course is ex- 
pected to be minimal with the periodic boundary conditions 
employed in the calculations. Therefore we expect the neigh- 
boring atoms to have significant 2p orbital contributions for 
concentrations well below ~ 1%. At higher concentrations, 
the contribution of neighboring atoms becomes less and less 
significant. A comparison of MPG and GA density of states in 
Fig.|5]indicates that the role of averaging is to slightly increase 
the contribution of p-orbitals which is essentially a statistical 
effect. 
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DISCUSSION AND CONCLUSION 

Examining the bond lengths show that Cu-C bond is elon- 
gated by approximately 0.2 relative to C-C bonds in undoped 
structure which has been almost compensated by a similar de- 
crease in the surrounding C-C bond lengths. C-C bonds far 
away from Cu undergo negligible changes of the order 0.01 A. 
Therefore in terms of structural changes, Cu substitution will 
give rise to local effects. In terms of electronic properties, in 
the regime of ^ 5% Cu substitution, the most significant pro- 
cess behind the formation of impurity band, is the hopping 
between Cu atoms themselves. In this regime a minimal two- 
band effective model can describe the electronic states near 
the Fermi energy; 

H ^ tddY^ {d\„d,„ + h.c.) - tpp ^ {p\aPj<y + h.c 

+tpd [d^Pja + h.c.) , 

where the hopping matrix elements tdd, tpp > are between 
the t2g hole bands and Ap electron bands of the adjacent Cu 
atoms residing at random sites If the inter-band ma- 

trix element tpd is negligible, one would have a quarter filled 
band of 3c? holes, along with a half-filled 4p band constitut- 
ing the impurity band. In this limit for large enough Coulomb 
correlations among the Wannier orbitals derived from parent 
4p states, the electron-like band undergoes a Mott transition, 
while the 1/4-filled hole bands derived from t2g states can be 
subject to various forms of density waves ll28ll . Theoretically 
speaking, in the opposite limit when tpd is large, the origi- 
nal half-filled electron band and the quarter filled hole band, 
get mixed giving rise to an effective quarter filled band whose 
electrons have a mixed character. 



(2) 
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This can be thought of as analogue of Zhang-Rice singlet 
state [230. Although this limit seems to be hard to realize 
in the experiment, but from theoretical point of view it may 
prove useful in studying transition from two-band situation to 
effective one-band model. In the intermediate regime where 
\tpd\ ^ tppjtdd^ copper doped diamond will offer an op- 
portunity to study both theoretically and experimentally the 
competition between the instabilities of a half-filled electronic 
band [29] with that of a quarter-filled hole band 12811 . Our 
model Hamiltonian ^ defines the competition between in- 
stabilities of two bands in 3D. The randomness in the matrix 
elements can not be avoided, as there is no control over the lo- 
cation of Cu impurities. Robustness of Lowdin charges in var- 
ious geometries imposed the restriction tdd, tpp > for study- 
ing Anderson localization effects in such 3D bands. The fact 
that in tpd ^ tpp, tdd limit electronic bands of Ap orbitals are 
half-filled, suggests the possibility of Mott- Anderson physics 
in Ap orbitals. 



To summarize, from both experimental and theoretical 
point of view, impurity band formation in Cu/C system pro- 
vides interesting ground to study the orbital selective phenom- 
ena and novel forms of Anderson localization properties in a 
combination of electron and hole like bands with different fill- 
ings in 3D. 



ACKNOWLEDGEMENT 

We thank Dr M. Alaei for fruitful discussions. H.H.F. 
would likes to thank R. Sarparast, K. Haghighi and M. Hafez 
for assistance in computing. S.A.J, was supported by the Vice 
Chancellor for Research Affairs of the Isfahan University of 
Technology (lUT) and the National Elite Foundation (NEF) 
of Iran. 



[1] Field, J. E., The Properties of Natural and Synthetic Diamond , 

Academic Press, 1992. 
[2] Davies, G. (Ed.) Properties and Growth of Diamond, INSPEC, 

1994. 

[3] Nebel, C. E., and J. Ristein, Thin-Film Diamond I, Elsevier, 1st 
Edition, 2003. 

[4] Lawrence S. Pan, Don R. Kania Diamond electronic properties 

and applications. Springer, 1994. 
[5] E. A. Ekimov, V. A. Sidorov, E. D. Bauer, N. N. Melnik, N. J. 

Curro, J. D. Thompson, and S. M. Stishov: Nature (London) 

428 (2004) 542. 

[6] Y. Takano, M. Nagao, I. Sakaguchi, M. Tachiki, T. Hatano, K. 

Kobayashi, H. Umezawa, and H. Kawarada: Appl. Phys. Lett. 

85 (2004) 2851. 
[7] E. Bustarret et. al. Nature (London) 444 (2006) 465. 
[8] X. Blase, C. Adessi, D. Connetable, Phys. Rev. Lett. 93 (2004) 

237004. 

[9] E. Bourgeois, X. Blase, App. Phys. Lett. 90 (2007) 142511. 
[10] R. Kalish in: The Physics of Diamond, lOS Press, 1997, Edited 

by A. Paoletti, and A. Tucciarone. 
[11] A. Deneuville in: The Physics of Diamond, lOS Press, 1997, 

Edited by A. Paoletti, and A. Tucciarone. 
[12] Dan Wu, et al, Phys. Rev. B 73 (2006) 012501. 
[13] G. Baskaran, Sci. Technol. Adv. Mat. 7 (2007) S49-S53. 
[14] K.-W. Lee, and W. E. Pickett, Phys. Rev. Lett. 93 (2004) 

237003; Phys. Rev. B 73 (2006) 075105. 
[15] Yu. G. Pogorelov, V. M. Loktev, Phys. Rev. B 72 (2005) 

075213. 

[16] J. Nakamura, et. al, Phys. Rev. B 70 (2004) 245111. 

[17] H. Mukuda et al, Phys. Rev. B 75 (2007) 033301. 

[18] T. Shirakawa, S. Horiuchi, Y. Ohta, H. Fukuyama, J. Phys. Soc. 

Jpn. 76 (2007)014711. 
[19] A.T. Collins in: The Physics of Diamond, lOS Press, 1997, 

Edited by A. Paoletti, and A. Tucciarone. 
[20] M. Alaei, S. A. Jafari, and H. Akbarzadeh, J. Phys. Chem. 

Solids, 69 (2008) 3283. 
[21] T. Ghodselahi, et. al, to be published 

[22] P.W. Anderson, The theory of .superconductivity in the high-Tc 

cuprates, Princeton University Press, 1997. 
[23] R C. Zhang, T. M. Rice, Phys. Rev. B 37 (1988) 3759. 



6 



[24] P. Giannozzi, et. al, J. Phys.: Condens. Matter 21 (2009) 
395502. 

[25] J. P Perdew et. al. Phys. Rev. Lett. 77 (1996) 3865. 
[26] R. Fletcher, Practical Methods of Optimization, Wiley, New 
York, 1987. 

[27] T. Oguchi : Science and Technology of Advanced Materials 7 



(2006) S67-S70 
[28] H. Fukuyama, J. Phys. Soc. Jpn. 75 (2006) 051001. 
[29] M. Imada, A. Fujimori, Y. Tokura, Rev. Mod. Phys. 70 (1998) 

1039. 



